Valence-Bond-Solid state entanglement in a 2-D Cayley tree 



Heng Fan 1 , Vladimir Korepin 2 , Vwani Roychowdhury 1 
1 Electrical Engineering Department, University of California at Los Angeles, Los Angeles, CA 90095, USA 
2 C.N. Yang Lnstitute for Theoretical Physics, State University of New York at Stony Brook, Stony Brook, NY 11794-3840, USA 

(Dated: February 1, 2008) 

The Valence-Bond-Solid (VBS) states are in general ground states for certain gapped models. 
We consider the entanglement of VBS states on a two-dimensional Cayley tree. We show that the 
entropy of the reduced density operator does not depend on the whole size of the Cayley tree. We 
also show that asymptotically, the entropy is liearly proportional to the number of singlet states cut 
by the reduced density operator of the VBS state. 
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PACS numbers: 75.10.Pq, 03.67.Mn, 03.65.Ud, 03.67.-a 



The entanglement of quantum states, and in particu- 
lar of the ground states, related with spin systems has 
been attracting a great deal of interest, see for example 

[aiiaiiiiaiiiiiaiiiiiiiiiiiSiiiiiiiiiia. 

We can quantify the entanglement by the von Neumann 
entropy of the reduced density operator of the ground 
states. This quantity for discrete spin chain and lattice 
models is analogous to the geo metric entropy in the con- 
tinous field theory [l^. l2fl l2l|. The geometric entropy 
plays an important role in the quantum field theory and is 
considered to be related to the Bekenstein-Hawking black 
hole entropy |2^| . It has been suggested flfl lifl ] that the 
geometric entropy is proportional to size of the boundary 
of the block. This problem has been recently studied for 
discrete case in the real free Klein- Gordon fields [lOj, and 
again, the entropy has been found to be related to the 
surface area of the fields. 

The entropy of a block of contiguous spins has been 
obtained for various one-dimensional spin chains. The 
higher-dimensional spin-chain case is in general more 
complicated, and only recently the entropy of the reduced 
density operator of ground states for higher-dimensional 
systmes was studied in 0, 0, 0, 0, [r| . For discrete 
2-D cases, the entropy is found to be linearly propor- 
tional to the boundary size of the reduced density oper- 
ator in the lattice|lfllll|. The fermion case is also stud- 
ied in [I2I lliij . In this Letter, we consider the entropy 
of a 2-D Valence-Bond Solid (VBS) state on a Cayley 
tree, which is the ground state of the Ameck-Kenncdy- 
Lieb-Tasaki (AKLT) model^H- The AKLT model is 
a gapped model [23 and has been well-studied in con- 
densed matter physics. Certain entanglement properties 
of the ground state of the AKLT model were already 
studiedlll But the entropy of a block of 

spins in a 2-D Cayley tree case is not available. There 
are several reasons that make this quantity of interest 
to the condensed matter community. First, the entropy 
quantifies the entanglement between the spins in the re- 
duced density operator with the rest of the ground state. 
This is important for qunatum computation and quan- 
tum information, and VBS states for quantum computa- 
tion are presnted in Ref . |2fij . Secondly, it is interest- 



ing to know whether this quantity is related to some 
macroscopic properties, such as susceptibility as pointed 
out in Ref.|27|. Third, it is generally expected that the 
entanglement inherent in a system is possibly responsi- 
ble for the phenomena of quantum phase-transitions 0. 
Since the AKLT model is a well studied model in con- 
densed matter physics, it is interesting to uncover rela- 
tionships between the entropy of a block of spins in the 
ground state and other measurable physical quantities, 
such as the correlation functions. Finally, the entangle- 
ment properties can indicate whether the density matrix 
renormalization method |28j can be used to efficiently 
simulate the quantum many-body svstems [2^.l3fj| . 
The Hamiltonian of the AKLT model is written as 231 
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where Sj is the spin operator on lattice site i, P(S) is the 
the orthogonal projection, and are undordered pairs 
in the lattice. The ground states of the AKLT model are 
known as VBS states, which are constructed from the 
singlet state \^f~) = 4=(| fj.) — | It). For convenience 

sake, we replace the singlet state by |W) = tt) + I 11 
). This substituition does not change any properties of 
the entanglement of the VBS state in this paper. 

Let's consider a C4 Cayley tree as in Ref.|23j and as 
showed in the Figure 1. At each site, three singlet states 
are conneted to it. Each dot represents a spin-i. A 
symmetrization of the three spin-l/2's creates a spin-3/2 
at each site. First, we only onsider one lattice site A, and 
three spin 1/2's (Ai,A2, A3) are located at this site. The 
three singlet states are \^!)a 1 b 1 \ 1 &)a 2 c 2 \^)a 3 d 3 ', see Fig 
2. 

For the VBS state, we need to consider a symmetriza- 
tion, Pa, at each lattice site A. We realize the sym- 
metrization operator Pa as follows: 



Pa 



i(|3 T>(3 T 



2U>(2U| 
+1 T,2|)(T,2|| + |3|)(3| |), (2) 
where state \i J.) is a symmetrized state with i spin 



up, and j spin down. For example |2 f, |) 
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FIG. 1: Cayley tree: each lattice site has three singlet states 
conneted to it, and a symmetrization operator is performed 
on each lattice site. 




FIG. 2: One site reduced density operator. Site A is four 
dimension which is a symmetrization of three spin- 1/2 at 
Ai,A2, A3. Site A is one part of three singlet states. 

) + I IT 4) + I lit))- For convenience, we also will use 
the notations |3 T) = |4), |2 U) = |3), | |,2 |) = |2), 
|3i) = |l). 

One site reduced density operator: Let's start from 
lattice site A, three singlet states are conneted with site 
A, \^) Ai b 1 \*)a 2 c 2 \^)a 3 d 3 , see Figure 2. If we choose 
symmetrizatin on A, we can find the following 

\®)a 123 b iC2 d 3 = Pa\*)a iBi \*)a 2 c 2 \*)a 3 d 3 

= ±(|44) + |33) + |22) + |11)), (3) 

where one site is lattice site A which is a symmetriza- 
tion on spin-1/2 sites ^1,^2,^3, other sites are bound- 
ary sites Bi,C2,D 3 . We can verify that even if we just 
project the state at lattice site A to the symmetric sub- 
space, then automatically, the rest of the sites B\, C2, D3 
are also projected to the corresponding symmetric sub- 
spaces. The reduced density operator of lattice site A 
can be found easily: pa — jl, which is the identity oper- 
ator, /, in SU (4) with a normalization factor so that the 
trace is 1. For the Cayley tree case presented in Figure 
1 , the reduced density operator on each lattice site is the 
same except for the boundary sites. The von Neumann 
entropy of the reduced density operator is S{pa) = 2. 
So, we know that for the Cayley tree case, each lattice 
site is maximally entangled with the rest of the sites, and 
the entanglement is 2 ebits. This result can be directly 
extended to other types of Cayley tree cases, and we can 
show that each lattice site is maximally entangled with 
the rest of the lattice sites. 

The entropy of the reduced density operators does not 
depend on the whole size of the Cayley tree: Let's study 
the quantum state of ©■ This is a six-partite state, and 
we consider the entanglement across B\ : Ai^fiCiDa cut. 
It follows directly from the state in © that the state 



across the cut B\ : A\, 2,36*2-03 cut is maximally entan- 
gled. So, we can treat the quantum state | , F)bi: j 4i23C2-D3 
in © as a singlet state as follows, 

\*)b i: a 123 c 2 d 3 = (I T *T>fl 1 B 1 + I I *;) Bl sJ/V2, (4) 
where we denote 

l*T>a = ~j=\ TT)c 2 a j |4)a + -^(| U)c 2 d 3 

+ 1 11)c 2 d 3 )\S)a + ^\ H)c 2 d 3 \2)a, 

1*1)3! = 4=1 TT)c 2 u 3 |3)^ + -^=(| U)c 2 d 3 

+ 1 U)c 2 d 3 )\2)a + ^\ U)c 2 d 3 \1)a. (5) 

This can also be understood as the following: First, we 
have a singlet state shared by B1A1. To expand the Cay- 
ley tree from A\ to A123C23 (i.e., one leg is expanded to 
two legs) , we put two additional singlet states shared by 
A 2 , C2 and A 3 , C3, and perform the symmetrical projec- 
tion on site A 12 3- The final result is that we just replace 
the spin up and spin down in A% by \^>f) and Sim- 
ilarly, we can further expand the Cayley tree from sites 
C2, and C3. If we only consider the bipartite entangle- 
ment of the VBS state on a Cayley tree, the bipartite 
state B\B\ is just a singlet state no matter how many 
legs are represented by B\. 

In this Letter, we will consider the entropy of re- 
duced density operators of the VBS on a Cayley tree. 
The reduced density operators are one-site, 4-site,..., 
in size. In Figure 1, the cuts are shown 
as dashed circles. The entropy of the reduced density op- 
erator is the bipartite entanglement across these cuts, i.e., 
between the spins in the reduced density operator and the 
rest. 

Based on our observations, we present one of our main 
conclusions in this Letter: The entropy of reduced density 
operator S(pi) does not depend on the whole size of the 
Cayley tree, where i — 1,4, .... The reason follows from 
our observation in Eq.Q, which implies that for the re- 
duced density operator pi , the legs connected to it in the 
Cayley tree can be dealt as just singlet states, no matter 
how large the Cayley tree itself is. We can repeatedly use 
this relation and expand the Cayley tree, but the quan- 
tity S(pi) does not change. So, to study the quantity 
S(pi), we only need to consider the smallest Cayley tree 
since the entropy does not denpend on the size of the 
Cayley tree. For example, if we consider S(p4,), we will 
just need to study the six legs case as in Figure 3. This 
result is similar to the result for 1-D case0| where a dif- 
ferent proof is given. We should note that the restriction 
on i can be relaxed, even though we only consider special 
cases in this Letter. 

Entropy of 4-site reduced density operator: For a bi- 
parite pure state \^ab), we know S(pa) = S(pb), where 
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FIG. 3: 4 sites Cay ley tree. To study the entropy of 4-site 
reduced density operator, we can choose the smallest Cayley 
tree. There are 6 legs (singlet states) connected to the 4-site 
reduced density operator 




FIG. 4: Entropy of reduced density operator of two legs 
(dashed circle) is log 3 

Pa(b) is the reduced density operator. To calculate the 
entropy of the 4-site reduced density operator, we can 
calculate the entropy of the reduced density operator of 
six legs; see Figure 3. First we have the state a^c-iD?) 
which takes the form as in Figure 2. Then we can sub- 
stitute the spin up and spin down in Bi,C 2 ,D 3 by |\I>j) 
and and now the correponding Cayley tree reduces 
to that shown in Figure 3. After tracing over spin-2/3 
sites in ABCD, we obtain the reduced density operator 
as follows: 

P4 = V 3 + |pf 3 + ~pf 3 + X, (6) 

where I = i(pf + Pi) = diag. (1/3, 1/3, 1/3), p\ = 
diofl.(l/2,l/3,l/6), Pi = d»a5.(l/6,l/3,l/2), the basis 
chosen here is{|2t),|t,j),|2|)} which are basis vectors 
of the symmetrical subspace. Moreover, p-\(p{) is derived 
from | v I / t)(|4'j,)) by taking the trace on the spin 2/3 site, 

X = -!- [I ® (x ® x l + x* <g> a:) + (x ® x t + x* ® x) ® I 
108 

+x<g)I(g>x t + x t ®I®x], (7) 

where x = |2 TXT, 4 I + I T4)(2 I |. Since the entropy 
of the density operator of ABCD is equal to six legs 
conneted to it, S(p4,) = S (pi), the entropy of the 4-site 
reduced density operator is 

53 13 65 
5(P4)= 54 + T l0g3 ~ 108 1Og5w 4 - 735 ' (8) 

where log has the base 2. We can find that P4 is actu- 
ally like I® 3 which has the entropy 3 log 3 « 4.755. The 
difference between these two entropies is less than 0.5%. 
In fact, the off-diagonal entries X in I© are small, and 
2 J2 a =u Pa 3 ^ S c l° se to ^® 3 - For simplicity, we can use 
the fidelity to define the distance between pj-^ and /. 
The fidelity is defined as F(pa,Pb) = Try/ ^Jp~APB\fpA- 



The fidelity is 0.977 for both pj and p^, and we know 
that pi 1 is like the identity /. That's the reason why 
S(p 4 )« / :i . 

General entropy of reduced density operator of the VBS 
state in Cayley tree: In the 2-D case presented in 
Ref. 0, 0| , the entropy of the reduced density operator 
is linearly proportional to size of the boundary. A toy 
model can be considered for VBS states: For the AKLT 
model, the ground states are constructed from the sin- 
glet states by some projection operators. Suppose we 
consider a state without any projections, the entropy of 
this state is the number of singlet states that cross the 
boundary. The result that the entropy is linearly propor- 
tional to the boundary size can be roughly understood as 
that the entropy is linearly proportional to the number of 
states (like the singlet states) cut by the boundary. For 
the VBS state in a Cayley tree, we expect similar results. 
In the last section, we found that S(pi) » 3 log 3. The 
number of singlet states crossing the boundary is 6. We 
can expect that in general we have S(pi) ss ^log^, i.e., 
N singlet states cross the boundary (i — ^ fe 3 x 2 k + 1, 
here N and i has a restriction between them). 

Let's show next the result: The entropy of the reduced 
density operator S(pi) is upper bounded by ^logS. A 
simple observation is that the upper bound is iV since N 
singlet states cross the boundary and the projection can 
only reduce this quantity. We next show that a tighter 
upper bound can be found. Consider the state presented 
in Q, the reduced density operator crossing two singlet 
states is diag. (1/3, 1/3, 1/3), and the entropy is log3, see 
Figure 4. This can be considered as a building block for 
the Cayley tree, thus we know that S(pi) < ^ log 3. 

This is just the first order upper bound, and we can 
further consider a second order upper bound: To find the 
entropy of 4 legs, for this case, we suppose the Cayley 
tree is relatively large. Before we present this second 
order upper bound, we present the techniques used in 
our calculations by the following example. 

We use pio as an example. Starting from p± in © , we 
can expand the Cayley tree by substituting the state | f ) 
and I I) by |$"t) and l^i), respectively. By tracing out 
the spin-2/3 sites, we can obtain p 10 from p 4 in JfJJ by 
the following substitutions, 

p a - \l m + \pT+Y, a=T,i, (9) 

x -> j^(x® I + I<8> x), (10) 

Y = -^-(x ® x* + x* ® x). (11) 
108 

The final result shows that pio ~ I® 6 , and thus S(pia) w 
61og3. Actually, we find S(p 10 ) » 9.4891 » 6 log 3(1- e), 
and e is small and is about 0.22%. The fidelity between 
the term related with p a in pio with identity is 0.994 
which is better than the original 0.977 in p 4 case. The 
Eas. H9ll0lll|) and Eq.© provide us with an algorithm 
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to find the general reduced density operators of the VBS 
state on a Cayley tree. 

For the general case, we can expand the Cayley tree 
and obtain the reduced density operator pi by relations 
Pl()lll|l . Suppose i = J2tLo 3 x 2 k + 1; then the general 
form will be as follows: 

ft = 2 -I^ M - 1 + \ £ (f(p a ,I))+g(x,I) 

» J®3x2 M -\ (12) 

The last equation follows from the fact that the second 
term is more like the identity operator as i becomes large, 
and the third term remains small and can be omitted. We 
know that S(pi) « 3 x 2 M - x log3,M > 1. The reduced 
density operator cuts N = 3 x 2 M singlet states. And 
we thus estimate that asymptotically the entropy of VBS 
in a Cayley tree is linearly proportional to the number of 
singlet states across the boundary. 

Wee should point out that though S(pi) ~ 3 x 
2 M_1 log3,M > 1, the upper bound 3 x 2 M " 1 log3 can 
never be saturated. The reason is that, as we pointed out, 
this is the first order upper bound. We can also compute 
a second order upper bound. Using the same method as 
presented above, we can obtain the reduced density oper- 
tor of 4 legs: Start from the 2-leg reduced density opera- 
tor |2 t>(2 T | + | T4X1U | + |21)<2 1 |, replace the state 
| t) by ^j) and similarly for | J.), trace out the corre- 
sponding sites, and then we can find the entropy of 4 legs 
to be S4 « 3.1631. We know the first order upper bound 
should give 2 log 3 w 3.1699. The difference between 
these two bounds is about 0.21 percent. So we have a 
more exact result: S(pi) < 3 x 2 M ~ 2 x 3.1631. Of course, 
we can also compute higher-order upper bounds. But as 
in 1-D VBS state in [Pja], the correction of higher order 
bounds will decay exponentially. We performed the cal- 
culations for the density operator of 8 legs and 16 legs cor- 
responding to third and fourth order corrections for the 
upper bound. As expected, we get 58/4 log 3 = (1 — e') 
and S'i6/81og3 = (1 — e"), both e' and e" are around 
0.21% which have almost no difference with the second 
order correction. So roughly, the first order upper bound 
log 3 still works. 

Observing that the result of S(pio) actually provides 
a lower-bound correction for the general S(pi), when i is 
large, asympototically we have 

1 - ei < S{ Pi )/3 x 2 M -* log 3 < 1 - e u (13) 

where safely we can set e; = 0.23%, and e u = 0.20%. 
Now we summarize our main result: Asymptotically, the 
entropy of the VBS state in a Cayley tree is linearly pro- 
portially to the number of singlet states that cross the 
boundary. We estimate that asymptotically e « 0.22%. 

We estimated the entropy of the reduced density oper- 
ators that are circles, as presented in Figure 1. But our 



method works for all kinds of reduced density operators 
on the Cayley tree. The only difference is that the sub- 
stitutions I|9I10J1 will depend on the form of the reduced 
density operators in the Cayley tree. We expect that our 
result that asymptotically, the entropy is liearly proporti- 
nal to the number of singlet states across the boundary, 
still holds given that the number of cut singlet states is 
large. We performed the calculations for a density oper- 
ator with 16 legs (not a circle) extended from p w , and 
we found that the Eq. (|13|l still holds for this case. Some 
other forms, with fewer legs, have been also checked and 
our conjecture holds for all these cases as well. 

In fact, we actually provide an algorithm to find the 
reduced density operators of VBS state on a Cayley tree. 
It will be interesting to apply this algorithm in the sim- 
ulations of quantum many-body systems. 
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